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ABSTRACT 

Turbulence is the most common state of astrophysical flows. In typical astrophysical 
fluids, turbulence is accompanied by strong magnetic fields, which has a large impact 
on the dynamics of the turbulent cascade. Recently, there has been a significant break¬ 
through on the theory of magnetohydrodynamic (MHD) turbulence. For the first time 
we have a scaling model that is supported by both observations and numerical simu¬ 
lations. We review recent progress in studies of both incompressible and compressible 
turbulence. We compare Iroshnikov-Kraichnan and Goldreich-Sridhar models, and dis¬ 
cuss scalings of Alfven, slow, and fast waves. We discuss the implications of this new 
insight into MHD turbulence for cosmic ray transport. 


1. Introduction 

Most astrophysical systems, e.g. accretion disks, stellar winds, the interstellar medium (ISM) 
and intercluster medium are turbulent with an embedded magnetic field that influences almost all of 
their properties. This turbulence which spans from km to many kpc (see discussion in Armstrong, 
Rickett, Spangler 1995; Scalo 1987; Lazarian, Pogosyan, & Esquivel 2002) holds the key to many 
astrophysical processes. For instance, propagation of cosmic rays and their acceleration is strongly 
affected by MHD turbulence. Recent research has shown that a substantial part of the earlier 
results in the field require revision. Earlier research used ad hoc models of MHD turbulence and 
this entailed erroneous conclusions. 

Before we start a discussion of MHD turbulence let us recall some basic properties of the 
hydrodynamic turbulence. All turbulent systems have one thing in common: they have a large 
“Reynolds number” (Re = LV/v ; L= the characteristic scale or driving scale of the system, V=the 
velocity difference over this scale, and i/=viscosity), the ratio of the viscous drag time on the 
largest scales (L 2 /v) to the eddy turnover time of a parcel of gas (L/V). A similar parameter, the 
“magnetic Reynolds number”, Rm (= LV/ry r/=magnetic diffusion), is the ratio of the magnetic 
field decay time (L 2 jrf) to the eddy turnover time (L/V). The properties of the flows on all scales 
depend on Re and Rm. Flows with Re < 100 are laminar; chaotic structures develop gradually 
as Re increases, and those with Re 10 3 are appreciably less chaotic than those with Re 10 7 . 
Observed features such as star forming clouds and accretion disks are very chaotic with Re > 10 8 
and Rm > 10 16 . 
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Let us start by considering incompressible hydrodynamic turbulence, which can be described 
by the Kolmogorov theory (Kolmogorov 1941). Suppose that we excite fluid motions at a scale 
L. We call this scale the energy injection scale or the largest energy containing eddy scale. For 
instance, an obstacle in a flow excites motions on scales of the order of its size. Then the energy 
injected at the scale L cascades to progressively smaller and smaller scales at the eddy turnover rate, 
i.e. t^ 1 & vi/l, with negligible energy losses along the cascade 1 . Ultimately, the energy reaches 
the molecular dissipation scale Id, i.e. the scale where the local Re ~ 1, and is dissipated there. 
The scales between L and Id are called the inertial range and it typically covers many decades. 
The motions over the inertial range are self-similar and this provides tremendous advantages for 
theoretical description. 

The beauty of the Kolmogorov theory is that it does provide a simple scaling for hydrodynamic 
motions. If the velocity at a scale l from the inertial range is vi, the Kolmogorov theory states that 
the kinetic energy (pvf ~ v? as the density is constant) is transferred to next scale within one eddy 
turnover time (l/vi). Thus within the Kolmogorov theory the energy transfer rate (vf /(l/vi)) is 
scale-independent, 

V ^ 

~ — j— = constant, (1) 

teas (l/Vl) 

and we get the famous Kolmogorov scaling 

vi oc Z 1/3 . (2) 

The one-dimensional 2 energy spectrum E(k ) is the amount of energy between the wavenumber 
k and k + dk divided by dk. When E(k) is a power law, kE{k) is the energy near the wavenumber 
k ~ l/l. Since v‘f ~ kE{k), Kolmogorov scaling implies 

E(k) oc k~ 5/3 . (3) 

Kolmogorov scalings were the first major advance in the theory of incompressible turbulence. 
They have led to numerous applications in different branches of science (see Monin &; Yaglom 
1975). However, astrophysical fluids are magnetized and the a dynamically important magnetic 
field should interfere with eddy motions. 

Paradoxically, astrophysical measurements are consistent with Kolmogorov spectra (see LPE02). 
For instance, interstellar scintillation observations indicate an electron density spectrum is very close 
to —5/3 for 10 8 cm - 10 15 cm. (see Armstrong et al. 1995). At larger scales LPE02 summarizes the 
evidence of -5/3 velocity power spectrum over pc-scales in HI. Solar-wind observations provide in- 
situ measurements of the power spectrum of magnetic fluctuations and Leamon et al. (1998) also 

lr This is easy to see as the motions at the scales of large eddies have Re 1. 

2 Dealing with observational data, e.g. in LPE02 (Lazarian et al. 2002), we deal with three dimensional energy 
spectrum P(k), which, for isotropic turbulence, is given by E{k) = 4ivk 2 P(k). 
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obtained a slope of ~ -5/3. Is this a coincidence? What properties is the magnetized compressible 
ISM expected to have? We will deal with these questions, and some related issues, below. 

Here we discuss a focused approach which aims at obtaining a clear understanding on the 
fundamental level, and considering physically relevant complications later. The creative synthesis 
of both approaches is the way, we think, that studies of astrophysical turbulence should proceed 3 . 
Certainly an understanding of MHD turbulence in the most ideal terms is a necessary precursor to 
understanding the complications posed by more realistic physics and numerical effects. For review 
of general properties of MHD, see a recent book by Biskanrp (1993). 

In what follows, we first consider observational data that motivate our study (§2), then discuss 
theoretical approaches to incompressible MHD turbulence (§3). We move to the effects of com¬ 
pressibility in §4 and discuss implications of our new understanding of MHD turbulence for cosmic 
ray dynamics in §5. We present the summary in §6. 


2. Observational Data 

Kolmogorov turbulence is the simplest possible model of turbulence. Since it is incompressible 
and not magnetized, it is completely specified by its velocity spectrum. If a passive scalar field, 
like “dye particles” or temperature inhonrogeneities, is subjected to Kolmogorov turbulence, the 
resulting spectrum of the passive scalar density is also Kolmogorov (see Lesieur 1990; Warhaft 2000). 
In compressible and magnetized turbulence this is no longer true, and a complete characterization 
of the turbulence requires not only a study of the velocity statistics but also the statistics of density 
and magnetic fluctuations. 

Direct studies of turbulence 4 have been done mostly for interstellar medium and for the Solar 
wind. While for the Solar wind in-situ measurements are possible, studies of interstellar turbulence 
require inverse techniques to interpret the observational data. 

Attempts to study interstellar turbulence with statistical tools date as far back as the 1950s 
(von Horner 1951; Kampe de Feriet 1955; Munch 1958; Wilson et al. 1959) and various directions 
of research achieved various degree of success (see reviews by Kaplan & Pickelner 1970; Dickman 
1985; Armstrong et al. 1995; Lazarian 1999a, 1999b; LPE02). 


3 Potentially our approach leads to an understanding of the relationship between motions at a given time at small 
scales (subgrid scales) and the state of the flow at a previous time at some larger, resolved, scale. This could lead to 
a parametrization of the subgrid scales and to large eddy simulations of MHD. 

4 Indirect studies include the line-velocity relationships (Larson 1981) where the integrated velocity profiles are 
interpreted as the consequence of turbulence. Such studies do not provide the statistics of turbulence and their 
interpretation is very model dependent. 
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2.1. Solar wind 

Solar wind (see review Goldstein & Roberts 1995) studies allow pointwise statistics to be 
measured directly using spacecrafts. These studies are the closest counterpart of laboratory mea¬ 
surements. 

The solar wind flows nearly radially away from the Sun, at up to about 700 km/s. This is much 
faster than both spacecraft motions and the Alfven speed. Therefore, the turbulence is “frozen” 
and the fluctuations at frequency / are directly related to fluctuations at the scale k in the direction 
of the wind, as k = 2irf /v, where v is the solar wind velocity (Horbury 1999). 

Usually two types of solar wind are distinguished, one being the fast wind which originates in 
coronal holes, and the slower bursty wind. Both of them show, however, / -5 / 3 scaling on small 
scales. The turbulence is strongly anisotropic (see Klein et al. 1993) with the ratio of power in 
motions perpendicular to the magnetic field to those parallel to the magnetic field being around 
30. The intermittency of the solar wind turbulence is very similar to the intermittency observed in 
hydrodynamic flows (Horbury & Balogh 1997). 


2.2. Electron density statistics 

Studies of turbulence statistics of ionized media (see Spangler & Gwinn 1990) have provided 
information on the statistics of plasma density at scales 10 8 -10 15 cm. This was based on a clear 
understanding of processes of scintillations and scattering achieved by theorists 5 (see Narayan &; 
Goodman 1989; Goodman & Narayan 1985). A peculiar feature of the measured spectrum (see 
Armstrong et al. 1995) is the absence of the slope change at the scale at which the viscosity by 
neutrals becomes important. 

Scintillation measurements are the most reliable data in the “big power law” plot in Armstrong 
et al. (1995). However there are intrinsic limitations to the scintillations technique due to the 
limited number of sampling directions, its relevance only to ionized gas at extremely small scales, 
and the impossibility of getting velocity (the most important!) statistics directly. Therefore with 
the data one faces the problem of distinguishing actual turbulence from static density structures. 
Moreover, the scintillation data does not provide the index of turbulence directly, but only shows 
that the data are consistent with Kolmogorov turbulence. Whether the (3D) index can be -4 instead 
of -11/3 is still a subject of intense debate (Higdon 1984; Narayan & Goodman 1989). In physical 
terms the former corresponds to the superposition of random shocks rather than eddies. 

Additional information on the electron density is contained in the Faraday rotation measures 
of extragalactic radio sources (see Simonetti & Cordes 1988; Simonetti 1992). However, there is 
so far no reliable way to disentangle contributions of the magnetic held and the density to the 


5 In fact, the theory of scintillations was developed first for the atmospheric applications. 
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signal. We feel that those measurements may give us the magnetic field statistics when we know 
the statistics of electron density better. 


2.3. Velocity and density statistics from spectral lines 

Spectral line data cubes are unique sources of information on interstellar turbulence. Doppler 
shifts due to supersonic motions contain information on the turbulent velocity field which is oth¬ 
erwise difficult to obtain. Moreover, the statistical samples are extremely rich and not limited 
to discrete directions. In addition, line emission allows us to study turbulence at large scales, 
comparable to the scales of star formation and energy injection. 

However, the problem of separating velocity and density fluctuations within HI data cubes 
is far from trivial (Lazarian 1995, 1999b; Lazarian & Pogosyan 2000; LPE02). The analytical 
description of the emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan 
(2000) (see also Lazarian 1999b; LPE02 for reviews) shows that the relative contribution of the 
density and velocity fluctuations depends on the thickness of the velocity slice. In particular, the 
power-law asymptote of the emissivity fluctuations changes when the dispersion of the velocity at 
the scale under study becomes of the order of the velocity slice thickness (the integrated width of the 
channel map). These results are the foundation of the Velocity-Channel Analysis (VCA) technique 
which provides velocity and density statistics using spectral line data cubes. The VCA has been 
successfully tested using data cubes obtained via compressible magnetohydrodynamic simulations 
and has been applied to Galactic and Small Magellanic Cloud atomic hydrogen (HI) data (Lazarian 
et al. 2001; Lazarian & Pogosyan 2000; Stanimirovic & Lazarian 2001; Deshpande, Dwarakanath, 
& Goss 2000). Furthermore, the inclusion of absorption effects (Lazarian & Pogosyan 2002) has 
increased the power of this technique. Finally, the VCA can be applied to different species (CO, 
H a etc.) which should further increase its utility in the future. 

Within the present discussion a number of results obtained with the VCA are important. First 
of all, the Small Magellanic Cloud (SMC) HI data exhibit a Kolmogorov-type spectrum for velocity 
and HI density from the smallest resolvable scale of 40 pc to the scale of the SMC itself, i.e. 4 kpc. 
Similar conclusions can be inferred from the Galactic data (Green 1993) for scales of dozens of 
parsecs, although the analysis has not been done systematically. Deshpande et al. (2000) studied 
absorption of HI on small scales toward Cas A and Cygnus A. Within the VCA their results can be 
interpreted as implying that on scales less than 1 pc the HI velocity is suppressed by ambipolar drag 
and the spectrum of density fluctuations is shallow P(k) ~ k~ 2 - 8 . Such a spectrum (Deshpande 
2000) can account for the small scale structure of HI observed in absorption. 
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2.4. Magnetic field statistics 

Magnetic field statistics are the most poorly constrained aspect of ISM turbulence. The po¬ 
larization of starlight and of the Far-Infrared Radiation (FIR) from aligned dust grains is affected 
by the ambient magnetic fields. Assuming that dust grains are always aligned with their longer 
axes perpendicular to magnetic field (see the review Lazarian 2000), one gets the 2D distribution 
of the magnetic field directions in the sky. Note that the alignment is a highly non-linear process 
in terms of the magnetic field and therefore the magnetic field strength is not available 6 . 

The statistics of starlight polarization (see Fosalba et al. 2002) is rather rich for the Galactic 
plane and it allows to establish the spectrum 7 E(K ) ~ K ~ 1 ' 5 , where K is a two dimensional wave 
vector describing the fluctuations over sky patch. 8 For uniformly sampled turbulence it follows 
from Lazarian & Shutenkov (1990) that E(K ) ~ K a for K < K 0 and K for K > Kq, where 
K- 1 is the critical angular size of fluctuations which is proportional to the ratio of the injection 
energy scale to the size of the turbulent system along the line of sight. For Kolmogorov turbulence 
a = -11/3. 

However, the real observations do not uniformly sample turbulence. Many more close stars are 
present compared to the distant ones. Thus the intermediate slops are expected. Indeed, Cho & 
Lazarian (2002b) showed through direct simulations that the slope obtained in Fosalba et al. (2002) 
is compatible with the underlying Kolmogorov turbulence. At the moment FIR polarimetry does 
not provide maps that are really suitable to study turbulence statistics. This should change soon 
when polarimetry becomes possible using the airborne SOFIA observatory. A better understanding 
of grain alignment (see Lazarian 2000) is required to interpret the molecular cloud magnetic data 
where some of the dust is known not to be aligned (see Lazarian, Goodman, & Myers 1997 and 
references therein). 

Another way to get magnetic field statistics is to use synchrotron emission. Both polarization 
and intensity data can be used. The angular correlation of polarization data (Baccigalupi et al. 
2001) shows the power-law spectrum A' -1 ' 8 and we believe that the interpretation of it is similar 
to that of starlight polarization. Indeed, Faraday depolarization limits the depth of the sampled 
region. The intensity fluctuations were studied in Lazarian & Shutenkov (1990) with rather poor 
initial data and the results were inconclusive. Cho & Lazarian (2002b) interpreted the fluctuations 
of synchrotron emissivity (Giardino et al. 2001, 2002) in terms of turbulence with Kolmogorov 
spectrum. 


®The exception to this may be the alignment of small grains which can be revealed by microwave and UV po¬ 
larimetry (Lazarian 2000). 

7 Earlier papers dealt with much poorer samples (see Kaplan & Pickelner 1970) and they did not reveal power-law 
spectra. 

8 This spectrum is obtained by Fosalba et al. (2002) in terms of the expansion over the spherical harmonic basis 
Yim.- For sufficiently small areas of the sky analyzed the multipole analysis results coincide with the Fourier analysis. 
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3. Theoretical Approaches to MHD Turbulence 

Attempts to describe magnetic turbulence statistics were made by Iroshnikov (1963) and 
Kraichnan (1965). Their model of turbulence (IK theory) is isotropic in spite of the presence 
of the magnetic field. 

For simplicity, let us suppose that a uniform external magnetic held (Bo) is present. In the 
incompressible limit, any magnetic perturbation propagates along the magnetic held line. Since 
wave packets are moving along the magnetic held line, there are two possible directions for prop¬ 
agation. If all the wave packets are moving in one direction, then they are stable to nonlinear 
order (Parker 1979). Therefore, in order to initiate turbulence, there must be opposite-traveling 
wave packets and the energy cascade occurs only when they collide. The IK theory starts from this 
observation, one of the consequences of which is the modification of the energy cascade timescale: 
teas ~ (L/V)(Va/^ 0> where Va = B^/^Airp is Alfven velocity of the mean held. Here, the IK 
theory assumes that opposite-traveling isotropic wave packets of similar size interact. From this 
and the scale-invariance of energy cascade rate, they obtained 

Iroshnikov-Kraichnan: E(k) oc k~ 3 ^ 2 . (4) 

However, the presence of the uniform magnetic component has non-trivial dynamical effects 
on the turbulence fluctuations. One obvious effect is that it is easy to mix held lines in directions 
perpendicular to the local mean magnetic held and much more difficult to bend them. The IK theory 
assumes isotropy of the energy cascade in Fourier space, an assumption which has attracted severe 
criticism (Montgomery & Turner 1981; Shebalin, Matthaeus, & Montgomery 1983; Montgomery & 
Matthaeus 1995; Sridhar & Goldreich 1994; Matthaeus et al. 1998). Mathematically, anisotropy 
manifests itself in the resonant conditions for 3-wave interactions: 

ki + k 2 = k 3 , (5) 

wi + (j 2 = w 3 , (6) 

where k’s are wavevectors and u/s are wave frequencies. The first condition is a statement of 
wave momentum conservation and the second is a statement of energy conservation. Alfven waves 
satisfy the dispersion relation: u ; = V)i|£:|||, where fey is the component of wavevector parallel to the 
background magnetic held. Since only opposite-traveling wave packets interact, ki and k 2 must 
have opposite signs. Then from equations (5) and (6), either k» i or k^ 2 must be equal to 0 and 
Atm 3 must be equal to the nonzero initial parallel wavenumber. That is, zero frequency modes are 
essential for energy transfer (Shebalin et al. 1983). Therefore, in the wavevector space, 3-wave 
interactions produce an energy cascade which is strictly perpendicular to the mean magnetic held. 
However, in real turbulence, equation (6) does not need to be satished exactly, but only to within 
an an error of order 5lo ~ 1 /t cas (Goldreich & Sridhar 1995). This implies that the energy cascade 
is not strictly perpendicular to Bo, although clearly very anisotropic. 

We assume throughout this discussion that the rms turbulent velocity at the energy injection 
scale is comparable to the Alfven speed of the mean held and consider only scales below the energy 



injection scale. This is called strong turbulence regime. Note that, as a consequence, the regime of 
Bq 3> 5b is not considered in this review. However, the regime of Bq <C 5b is still relevant to the 
strong turbulence regime because scales below the energy equipartition scale is expected to fall in 
the strong turbulence regime (Cho & Vishniac 2000a). 

An ingenious model very similar in its beauty and simplicity to the Kolmogorov model has 
been proposed by Goldreich & Sridhar (1995; hereinafter GS95) for incompressible strong MHD 
turbulence. They pointed out that motions perpendicular to the magnetic field lines mix them 
on a hydrodynamic time scale, i.e. at a rate t~} s ~ k±vi, where k± is the wavevector component 
perpendicular to the local mean magnetic field and l ~ fe _1 (~ kj 1 ). These mixing motions couple 
to the wave-like motions parallel to magnetic field giving a critical balance condition 

k\\V A ~ k±v k , (7) 

where fey is the component of the wavevector parallel to the local magnetic field. When the typical 
feii on a scale k± falls below this limit, the magnetic field tension is too weak to affect the dynamics 
and the turbulence evolves hydro dynamically, in the direction of increasing isotropy in phase space. 
This quickly raises the value of fey. In the opposite limit, when fen is large, the magnetic field tension 
dominates, the error 5lu in the matching conditions is reduced, and the nonlinear cascade is largely 
in the k± direction, which restores the critical balance. 

If conservation of energy in the turbulent cascade applies locally in phase space then the energy 
cascade rate (vf/t cas ) is constant: (vf)/(l/vi) = constant. Combining this with the critical balance 
condition we obtain an anisotropy that increases with decreasing scale 

fey oc k 2 / 3 , (8) 

and a Kolmogorov-like spectrum for perpendicular motions 

vi oc l 1 / 3 , or,E(k) oc k ± 5 ^ 3 , (9) 

which is not surprising since the magnetic field does not influence motions that do not bend it. 
At the same time, the scale-dependent anisotropy reflects the fact that it is more difficult for the 
weaker, smaller eddies to bend the magnetic field. 

GS95 shows the duality of motions in MHD turbulence. Those perpendicular to the mean 
magnetic field are essentially eddies, while those parallel to magnetic field are waves. The critical 
balance condition couples these two types of motions. 

Numerical simulations (Cho & Vishniac 2000b; Maron & Goldreich 2001; Cho, Lazarian, & 
Vishniac 2002) show reasonable agreements with the GS95 model. 


4. Compressible Turbulence 

For the rest of the review, we consider MHD turbulence of a single conducting fluid. While 
the GS95 model describes incompressible MHD turbulence well, no universally accepted theory 
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exists for compressible MHD turbulence despite various attempts (e.g., Higdon 1984). Earlier 
numerical simulations of compressible MHD turbulence covered a broad range of astrophysical 
problems, such as the decay of turbulence (e.g. Mac Low 1998; Stone, Ostriker, & Gamrnie 1998) or 
turbulent modeling of the ISM (see recent review Vazquez-Semadeni 2002; see also Passot, Pouquet, 
& Woodward 1988; Vazquez-Semadeni, Passot, & Pouquet 1995; Passot, Vazquez-Semadeni, &; 
Pouquet 1995; Vazquez-Semadeni, Passot, & Pouquet 1996 for earlier pioneering 2D simulations 
and Ostriker, Gammie, & Stone 1999; Ostriker, Stone, & Gamrnie 2001; Padoan et al. 2001; 
Klessen 2001; Boldyrev 2002 for recent 3D simulations). In what follows, we concentrate on the 
fundamental properties of compressible MHD. 


4.1. Alfven, slow, and fast modes 


Let us start by reviewing different MHD waves. In particular, we describe the Fourier space 
representation of these waves. The real space representation can be found in papers on modern 
shock-capturing MHD codes (e.g. Brio & Wu 1988; Ryu & Jones 1995). For the sake of simplicity, 
we consider an isothermal plasma. Figure 1 and Figure 2 give schematics of slow and fast waves. 
For slow and fast waves, Bo, Vk (oc £), and k are in the same plane. On the other hand, for Alfven 
waves, the velocity of the fluid element (vjJa is orthogonal to the Bq — k plane. 


As before, the Alfven speed is Va = Bq/ y/Airpo, where po is the average density. Fast and slow 
speeds are 


c f,s ~ 


1 

2 



— 4 a 2 Vj cos 2 0 



( 10 ) 


where 6 is the angle between Bo and k. See Table 1 for the definition of other variables. When 
/?(/? = Pgj Pb=2cl 2 /Vf; P g = gas pressure, Pb= magnetic pressure; hereinafter [3 = average /3 = 
Pg/Ps) goes to zero, we have 


c f « V A , 


Table 1: Notations for compressible turbulence 


Notation 

Meaning 

a, c s , cf, Va 

sound, slow, fast, and Alfven speed 

SV, (6V)„ ( SV) f , (JV)a 

random (mis) velocity 

Previously we used V for the rms velocity 

Vi, (vi)s, ( vi)f , (vi)a 

velocity at scale l 

Vk, (Vk)s, (v k )/, (v k ) A 

velocity vector at wavevector k 

B 0 (=k„), k±, k, 0 , ... 

unit vectors 

is, if 

displacement vectors 
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(b) 

k ± 

k ± 

(c) 


Fig. 1.— (a) Directions of fast and slow basis vectors. £/ and represent the directions of 
displacement of fast and slow modes, respectively. In the fast basis (£/•) is always between k and 
k^. In the slow basis (£ s ) lies between 8 and Bo- Here, 8 is perpendicular to k and parallel to 
the wave front. All vectors he in the same plane formed by Bo and k. On the other hand, the 
displacement vector for Alfven waves (not shown) is perpendicular to the plane, (b) Directions 
of basis vectors for a very small (3 drawn in the same plane as in (a). The fast bases are almost 
parallel to kj_. (c) Directions of basis vectors for a very high (3. The fast basis vectors are almost 
parallel to k. The slow waves become pseudo-Alfven waves. 


a cos 8. 


( 11 ) 


Figure 1 shows directions of displacement (or, directions of velocity) vectors for these three 
modes. We will call them the basis vectors for these modes. The Alfven basis is perpendicular to 
both k and Bo, and coincides with the azimuthal vector <f> in a spherical-polar coordinate system. 
Here hatted vectors are unit vectors. The fast basis lies between k and kj_: 




1-VD +P/2 


l + VD-P/2 |>||. 


k± 


1 2 


fc|ik|| + /cj_kj_, 


( 12 ) 


where D = (1 + /3/2) 2 — 2/3 cos 2 8 , and P is the averaged P (=P 9 /Pb)- The slow basis £ s lies between 
8 and Bq (=kii): 


is oc fciikii + 


fen 


l-VD-P/2 
i + Vd + P/2 Lfci. 


k±k ± . 


(13) 


The two vectors £/ and are mutually orthogonal. Proper normalizations are required for both 
bases to make them unit-length. 

When p goes to zero (i.e. the magnetically dominated regime), becomes parallel to k^ and 
£ s becomes parallel to Bq (Fig. lb). The sine of the angle between Bq and is {P/2) sin 8 cos 8. 
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Fig. 2.— Waves in real space. We show the directions of displacement vectors for a slow wave 
(left) and a fast wave (right). Note that £ s lies between 9 and Bo (= ky) and £/ between k and kj_. 
Again, 9 is perpendicular to k and parallel to the wave front. Note also that, for the fast wave, 
for example, density (inferred by the directions of the displacement vectors) becomes higher where 
field lines are closer, resulting in a strong restoring force, which is why fast waves are faster than 
slow waves. 

When (3 goes to infinity (i.e. gas pressure dominated regime) 9 , £/ becomes parallel to k and £ s 
becomes parallel to 9 (Fig. lc). This is the incompressible limit. In this limit, the slow mode is 
sometimes called the pseudo-Alfven mode (Goldreich & Sridhar 1995). 


4.2. Theoretical considerations 

Here we address the issue of mode coupling in a low f3 plasma. It is reasonable to suppose that 
in the limit where /3 1 turbulence for Mach numbers (M s = 5V/a) less than unity should be 

largely similar to the exactly incompressible regime. Thus, Lithwick & Goldreich (2001) conjectured 
that the GS95 relations are applicable to slow and Alfven modes with the fast modes decoupled. 
They also mentioned that this relation can carry on for low f3 plasmas. For (3 S> 1 and M s > 1, we 
are in the regime of hydrodynamic compressible turbulence for which no theory exists, as far as we 
know. 


9 In this section, we assume that external mean held is strong (i.e. Va > (SV)) but finite, so that f3 —> oo means 
the gas pressure P g —> oo. 
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In the diffuse interstellar medium /3 is typically less than unity. In addition, it is ~ 0.1 or 
less for molecular clouds. There are a few simple arguments suggesting that MHD theory can be 
formulated in the regime where the Alfven Mach number (= 5V/Va ) is less than unity, although 
this is not a universally accepted assumption. Alfven modes describe incompressible motions. 
Arguments in GS95 are suggestive that the coupling of Alfven to fast and slow modes will be weak. 
Consequently, we expect that in this regime the Alfven cascade should follow the GS95 scaling. 
Moreover the slow waves are likely to evolve passively (Lithwick Sz Goldreich 2001). For a <C Va 
their nonlinear evolution should be governed by Alfven modes so that we expect the GS95 scaling 
for them as well. The phase velocity of Alfven waves and slow waves depend on a factor of cos0 
and this enables modulation of the slow waves by the Alfven ones. However, fast waves do not 
have this factor and therefore cannot be modulated by the changes of the magnetic field direction 
associated with Alfven waves. The coupling between the modes is through the modulation of the 
local Alfven velocity and therefore is weak. 

For Alfven Mach number (Af4) larger than unity a shock-type regime is expected. However, 
generation of magnetic field by turbulence (Cho Sz Vishniac 2000a) is expected for such a regime. 
It will make the steady state turbulence approach Ma ~ l. 10 Therefore in Cho & Lazarian (2002a) 
we consider turbulence in the limit M s > 1, Ma < 1, and (5 < 1. For these simulations, we mostly 
used M s ~ 2.2, Ma ~ 0.7, and /3 ~ 0.2. The Alfven speed of the mean external field is similar to 
the rms velocity {Va = 1, SV ~ 0.7, a = \/0.1), and we used an isothermal equation of state. 

Although the scaling relations presented below are applicable to sub-Alfvenic turbulence, we 
cautiously speculate that small scales of super-Alfvenic turbulence might follow similar scalings. 
Boldyrev, Nordlund, & Padoan (2001) obtained energy spectra close to E(k ) ~ k~ 1 - 74 in solenoidally 
driven super-Alfvenic supersonic turbulence simulations. The spectra are close to the Kolmogorov 
spectrum (~ k~ 5 ^ 3 ), rather than shock-dominated spectrum (~ k~ 2 ). This result might imply 
that small scales of super-Alfvenic MHD turbulence can be described by our sub-Alfvenic model 
presented below, which predicts Kolmogorov-type spectra for Alfven and slow modes. 


4.3. Coupling of MHD modes and Scaling of Alfven modes 

Alfven modes are not susceptible to collisionless damping (see Spangler 1991; Minter & Span¬ 
gler 1997 and references therein) that damps slow and fast modes. Therefore, we mainly consider 
the transfer of energy from Alfven waves to compressible MHD waves (i.e. to the slow and fast 
modes). 

In Cho & Lazarian (2002a), we carry out simulations to check the strength of the mode¬ 
mode coupling. We first obtain a data cube from a driven compressible numerical simulation with 


10 We suspect that simulations that show super-Alfvenic turbulence is widely spread in the ISM might not evolve 
for a long enough time to reach the steady state. 
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Fig. 3.— (Left) Decay of Alfvenic turbulence. The generation of fast and slow waves is not efficient. 
ft ~ 0.2, M s ~ 3. (Right) The ratio of (5V) 2 to (5V)\. The stronger the external field (Bq) is, the 
more suppressed the coupling is. The ratio is not sensitive to ft. From Cho & Lazarian (2002a) 


Bo/^/4irpo = 1. Then, after turning off the driving force, we let the turbulence decay. We go 
through the following procedures before we let the turbulence decay. We first remove slow and fast 
modes in Fourier space and retain only Alfven modes. We also change the value of Bo preserving its 
original direction. We use the same constant initial density po for all simulations. We assign a new 
constant initial gas pressure P g 11 . After doing all these procedures, we let the turbulence decay. 
We repeat the above procedures for different values of Bq and P g . Fig. 3a shows the evolution of 
the kinetic energy of a simulation. The solid line represents the kinetic energy of Alfven modes. 
It is clear that Alfven waves are poorly coupled to the compressible modes, and do not generate 
them efficiently 12 Therefore, we expect that Alfven modes will follow the same scaling relation as 
in the incompressible case. Fig. 3b shows that the coupling gets weaker as Bq increases: 


(SV) 2 f (5 V) a 
(4)i “ B 0 


(14) 


The ratio of (5V) 2 to (5V)\ is proportional to (5V)\/B'q. 

This marginal coupling is in good agreement with a claim in GS95, as well as earlier numerical 
studies where the velocity was decomposed into a compressible component vp and a solenoidal 


11 The changes of both Bo and P g preserve the Alfven character of perturbations. In Fourier space, the mean 
magnetic Held (Bo) is the amplitude of k = 0 component. Alfven components in Fourier space are for k 7^ 0 and 
their directions are parallel/anti-parallel to 4 (= Bo x kx). The direction of 4 does not depend on the magnitude 
of Bo or P g . 

12 As correctly pointed out by Zweibel (this volume) there is always residual coupling between Alfven and com¬ 
pressible modes due to steepening of Alfven modes. However, this steepening happens on time-scales much longer 
than the cascading time-scale. 
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Fig. 4.— (a) Alfven spectra follow a Kolmogorov-like power law. (b) The second-order structure 
function (SF 2 = < v(x + r) — v(x) >) for Alfven velocity shows anisotropy similar to the GS95. 
Conturs represent eddy shapes. From Cho & Lazarian (2002a). 

component The compressible component is curl-free and parallel to the wave vector k in 
Fourier space. The solenoidal component is divergence-free and perpendicular to k. The ratio 
X = (8V)c/($V)s is an important parameter that determines the strength of any shock (Passot 
et al. 1988; Pouquet 1999). Porter, Woodward, & Pouquet (1998) performed a hydrodynamic 
simulation of decaying turbulence with an initial sonic Mach number of unity and found that 
X 2 evolves toward ~ 0.11. Matthaeus et al. (1996) carried out simulations of decaying weakly 
compressible MHD turbulence (Zank & Matthaeus 1993) and found that x 2 ~ O(Mf), where M s is 
the sonic Mach number. In Boldyrev et al. (2001) a weak generation of compressible components 
in solenoidally driven super-Alfvenic supersonic turbulence simulations was obtained. 

Fig. 4 shows that the spectrum and the anisotropy of Alfven waves in this limit are compatible 
with the GS95 model: 

Spectrum of Alfven Modes: E(k ) oc kf 5 ^ 3 , (15) 

2/3 

and scale-dependent anisotropy k» oc kf that is compatible with the GS95 theory. 

4.4. Scaling of the slow modes 

Slow waves are somewhat similar to pseudo-Alfven waves (in the incompressible limit). First, 
the directions of displacement (i.e. £ s ) of both waves are similar when anisotropy is present. The 
vector £ s is always between 9 and ky. In Figure 1, we can see that the angle between 9 and ky gets 
smaller when k\\ <C k±. Therefore, when there is anisotropy (i.e. k\\ <C k±), of a low (3 plasma 
becomes similar to that of a high /3 plasma. Second, the angular dependence in the dispersion 
relation c s ~ a cos 9 is identical to that of pseudo-Alfven waves (the only difference is that, in slow 
waves, the sound speed a is present instead of the Alfven speed Va)- 
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Fig. 5.— (a) Slow spectra also follow a Kolmogorov-like power law. (b) Slow modes show anisotropy 
similar to the GS95 theory. From Cho & Lazarian (2002a). 


Goldreich & Sridhar (1997) argued that the pseudo-Alfven waves are slaved to the shear-Alfven 
(i.e. ordinary Alfven) waves in the presence of a strong Bo, meaning that the energy cascade of 
pseudo-Alfven modes is primarily mediated by the shear-Alfveen waves. This is because pseudo- 
Alfven waves do not provide efficient shearing motions. Similar arguments are applicable to slow 
waves in a low /3 plasma (Cho & Lazarian 2002a) (see also Lithwick & Goldreich 2001 for high-/? 
plasmas). As a result, we conjecture that slow modes follow a scaling similar to the GS95 model 
(Cho & Lazarian 2002a): 

Spectrum of Slow Modes: E s (k) oc fc ± 5 ^ 3 . (16) 


Fig. 5a shows the spectra of slow modes. For velocity, the slope is close to —5/3. Fig. 5b shows 
the contours of equal second-order structure function (SF 2 ) of slow velocity, which are compatible 
with k\\ oc k 2 / 3 scaling. 

In low (3 plasmas, density fluctuations are dominated by slow waves (Cho & Lazarian 2002a). 
From the continuity equation p = pS7 ■ v 

upk = Pok • v fc , (17) 


we have, for slow modes, (pk) s ~ Po(^fc)s/ a - Hence, this simple argument implies 

(Se\ _ m. „ 


where we assume that (5V) S ~ {5V)a and M s is the Mach number. On the other hand, only a small 
amount of magnetic field is produced by the slow waves. Similarly, using the induction equation 
(whk = k x (Bo x Vfc)), we have 

= (19) 
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Fig. 6.— (a) The power spectrum of fast waves is compatible with the IK spectrum, (b) The mag¬ 
netic second-order structure function of fast modes shows isotropy.From Cho & Lazarian (2002a). 

which means that equipartition between kinetic and magnetic energy is not guaranteed in low 
(3 plasmas. In fact, in Fig. 5a, the power spectrum for density fluctuations has a much larger 
amplitude than the magnetic field power spectrum. Since density fluctuations are caused mostly 
by the slow waves and magnetic fluctuation is caused mostly by Alfven and fast modes, we do 
not expect a strong correlation between density and magnetic field, which agrees with the ISM 
simulations (Padoan & Nordlund 1999; Ostriker et al. 2001; Vazquez-Semadeni 2002). 


4.5. Scaling of the fast modes 

Figure 6 shows fast modes are isotropic. The resonance conditions for interacting fast waves 

are: 


U\-\-U>2 = a> 3 , (20) 

ki + k 2 = k 3 . (21) 


Since u oc k for the fast modes, the resonance conditions can be met only when all three k vectors 
are collinear. This means that the direction of energy cascade is radial in Fourier space, and we 
expect an isotropic distribution of energy in Fourier space. 


Using the constancy of energy cascade and uncertainty principle, we can derive an IK-like 
energy spectrum for fast waves. The constancy of cascade rate reads 

1 1- r.3~,2 

—— = -- = constant. (22) 

teas teas 

On the other hand, t cas can be estimated as 


Vk 


Vk 


(v • Vv) k Ep+ q =k kv v v q 


r\j 


(23) 
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If contributions are random, the denominator can be written by the square root of the number of 
interactions (VA 7) times strength of individual interactions (~ kvT) 13 • Here we assume locality 
of interactions: p ~ q ~ k. Due to the uncertainly principle, the number of interactions becomes 
M ~ k{Ak) 2 , where Ak is the typical transversal (i.e. not radial) separation between two wave 
vectors p and q (with p + q = k). Therefore, the denominator of equation (23) is (k(Ak) 2 ) 1 / 2 kv\. 
We obtain an independent expression for t cas from the uncertainty principle (t cas A lo ~ 1 with 
Auj ~ Ak(Ak/k)). From this and equation (23), we get t cas ~ tlis/(k 2 Vk), which yields 

teas ~ 1 /k A v 2 k . (24) 

Combining equations (22) and (24), we obtain v k ~ k 7 / 2 , or E? ( k ) ~ k 2 v^ ~ k 3 / 2 . This is 
very similar to acoustic turbulence, turbulence caused by interacting sound waves (Zakharov 1967; 
Zakharov & Sagdeev 1970; L’vov, L’vov, & Pomyalov 2000). Zakharov & Sagdeev (1970) found 
E(k) oc A: -3 / 2 . However, there is debate about the exact scaling of acoustic turbulence. Here we 
cautiously claim that our numerical results are compatible with the Zakharov & Sagdeev scaling: 

Spectrum of Fast Modes: E? (k) ~ k ~ 3 / 2 . (25) 


Magnetic field perturbations are mostly affected by fast modes (Cho &; Lazarian 2002a) when 
/3 is small: 


Fast: 


m f (tv)f 

(' SV) A ~ (SV)a’ 


(26) 


if {5V) A ~ (5V) S 


The turbulent cascade of fast modes is expected to be slow and in the absence of collisionless 
damping they are expected to propagate in turbulent media over distances considerably larger than 
Alfven or slow modes. This effect is difficult to observe in numerical simulations with A B ~ Bq. 
A modification of the spectrum in the presence of the collisionless damping is presented in Yan & 
Lazarian (2002). 


5. Implication for Cosmic Ray Propagation 

Many astrophysical problems require some knowledge of the scaling properties of turbulence. 
Therefore we expect a wide range of applications of the established scaling relations. Here we show 
how recent progress in understanding MHD turbulence affects cosmic ray propagation. 

The propagation of cosmic rays is mainly determined by their interactions with electromagnetic 
fluctuations in interstellar medium. The resonant interaction of cosmic ray particles with MHD 


13 To be exact, the strength of individual interactions is ~ kv\sui9, where 9 is the angle between k and Bo. Thus 
marginal anisotropy is expected. It will be investigated elsewhere. 
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turbulence has been repeatedly suggested as the main mechanism for scattering and isotropizing 
cosmic rays. In these analysis it is usually assumed that the turbulence is isotropic with a Kol¬ 
mogorov spectrum (see Schlickeiser & Miller 1998). How should these calculations be modified? 

Consider resonance interaction first. Particles moving with velocity v get into resonance with 
MHD perturbations propagating along the magnetic field if the resonant condition is fulfilled, 
namely, lo = kuvp+nfl (n = ±1, 2 ...), where ui is the wave frequency, H = flo/y is the gyrofrequency 
of relativistic particle, p = cos ct, where a is the pitch angle of particles. In other words, resonant 
interaction between a particle and the transverse electric field of a wave occurs when the Doppler 
shifted frequency of the wave in the particle’s guiding center rest frame ui gc = lo — k\\vp is a multiple 
of the particle gyrofrequency. 

For cosmic rays, h\vp 3> lo, so the slow variation of the magnetic field with time can be 
neglected. Thus the resonant condition is simply kuvp + nfl = 0. From this resonance condition, 
we know that the most important interaction occurs at fey = k res = Sl/vu. 

It is intuitively clear that resonant interaction of particles in isotropic and anisotropic turbu¬ 
lence should be different. Chandran (2001) obtained strong suppression of scattering by Alfvenic 
turbulence when he treated turbulence anisotropies in the spirit of Goldreich-Sridhar model of in¬ 
compressible turbulence. His treatment was improved in Yan Sz Lazarian (2002, henceforth YL02) 
who used a more rigorous description of magnetic field statistics. Moreover, they took into account 
CR scattering by compressible MHD modes and found that fast modes absolutely dominate cosmic 
ray scattering. In our description we shall follow YL02 treatment of the problem. 

We employ quasi-linear theory (QLT) to obtain our estimates. QLT has been proved to be a 
useful tool in spite of its intrinsic limitations (Chandran 2000; Schlickeiser & Miller 1998; Miller 
1997). For moderate energy cosmic rays, the corresponding resonant scales are much smaller than 
the injection scale. Therefore the fluctuation on the resonant scale 5B <C Ho even if they are 
comparable at the injection scale. QLT disregards diffusion of cosmic rays that follow wandering 
magnetic field lines (Jokipii 1966) and this diffusion should be accounted separately. Obtained by 
applying the QLT to the collisionless Boltzmann-Vlasov equation, the Fokker-Planck equation is 
generally used to describe the evolvement of the gyrophase-average distribution function /, 


dj_ 

dt 


±( D ^l + D 

" ' p,,, t is/ip 


dp 


dp 


s J_d_ 

+ 00 


dp J p 2 dp 


p 


D t + D d ± 

,ip dp +pp d P 


where p is particle momentum. The Fokker-Planck coefficients D pix , D^ ip , D pp are the fundamental 
physical parameter for measuring the stochastic interactions, which are determined by the electro¬ 
magnetic fluctuations (Schlickeiser & Achatz 1993). From Ohm’s Law E(k) = —(l/c)v(k) x B 0 , 
we can get the electromagnetic fluctuations from correlation tensors of magnetic and velocity fluc¬ 
tuations Cij, Kij. Here, 


< Bi{k)B*( k') > / Bq = <$(k - k')My(k) 
< Vi(k)B*(k') > /V a B 0 = <5(k - k!)C i:i (k) 
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< ^(kK(k') > fvl = 5(k - k')^'(k). 


(27) 


For Alfven modes, Clio, Lazarian and Vishniac (2002) obtained 

Kij( k) = CJijkT^expi-L^kJk 2 / 3 ), 


(28) 


where Iij = {5{j — kikj/k 2 ^} is a 2D matrix in x-y plane, fey is the wave vector along the local mean 
magnetic field, k± is the wave vector perpendicular to the magnetic field and the normalization 
constant C a = L^ 1 / 3 /67r. We assume that for the Alfven modes M- tJ = K tJ . Cij = ctM jj where the 
fractional helicity —1 < a < 1 is independent of k (Chandran 2000). According to Cho & Lazarian 
(2002a), fast modes are isotropic and have one dimensional spectrum E(k ) oc At 3 / 2 . In low f3 
medium, the velocity fluctuations are always perpendicular to Bo for all k, while the magnetic 
fluctuations are perpendicular to k. Thus Kjj, Mjj of fast modes are not equal, 


(29) 


where Jjj = kikj/k“j_ is also a 2D tensor in x — y plane 14 . In high (3 medium, the velocity 
fluctuations are radial, i.e., along the direction of k. Fast modes in this regime are essentially 
sound waves compressing magnetic field (GS95; Lithwick & Goldreich 2001, Cho & Lazarian, in 
preparation). The compression of magnetic field depends on plasma (3. The corresponding x-y 
components of the tensors are 


" Mjj (k) ' 

L " 1 / 2 

cos 2 0 

Cij{ k) 

Jijk-P 2 
8vr 3 

a cos 2 6 

Kij( k) 
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Mij( k) 
Cij{ k) 
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t-1/2 

-sin 2 6Jijk~ 

87T 


cos 2 0//3 
CJCOS0//3 1 / 2 

1 


(30) 


Adopting the approach in Schlickeiser & Achatz (1993), we can obtain the Fokker-Planck coefficients 
in the lowest order approximation of Va/c , 
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14 Apparently My, Cij are 3D matrixes. However, the third dimension is not needed for our calculations. Mij is 
different from that in Schlickeiser & Miller (1998). The fact that the fluctuations <5B in fast modes are in the k-B 
plane place another constrain on the tensor so that the term 8ij doesn’t exist. 
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(31) 


where k mtn = L x , k max = &o/vth corresponds to the dissipation scale, m = 7 mn is the relativistic 
mass of the proton, v± is the particle’s velocity component perpendicular to Bo, 0 = arctan^/A^), 
c,n = (x ± iy)/V 2 represent left and right hand polarization 15 . 


5.1. Scattering by Alfvenic turbulence 


Noticing that the integrand for small k± is substantially suppressed by the exponent in the 
anisotropic tensor (see Eq. (28)) so that the large scale contribution is not important, we can 
simply use the asymptotic form of Bessel function for large argument. Then if the pitch angle a is 
not close to 0, we can derive the analytical result for anisotropic turbulence (YL02), 
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(32) 


where L is the injection sale, k max = 52o /vt.h corresponds to the dissipation scale, T[a, z\ represents 
the incomplete gamma function. 


The scattering frequency v = 2H wt /(l — /x 2 ) is plotted for different models in Fig.(7a). It is 
clear that anisotropy suppresses scattering. Although our results are larger than those obtained in 
Chandran (2001) using an ad hoc tensor with a step function 16 , they are still much smaller than 
the estimates for isotropic model. Unless we consider very high energy CRs (> 10 8 GeV) with 
corresponding Larrnor radius comparable to the injection scale, we can neglect scattering by the 
Alfvenic turbulence. What is the alternative way to scatter cosmic rays? 


5.2. Scattering by fast modes 

For compressible modes we discuss two types of resonant interaction: gyroresonance and 
transit-time damping; the latter requires longitudinal motions. The contribution from slow modes 


15 For D^p, the expression is only true for Alfv\’en modes. There are additional compressional terms for compress- 
able modes. 

16 We counted only the resonant term in Chandran (2001). The nonresonant term is spurious as noted in Chandran 
( 2001 ) 
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is not larger than that by Alfven modes since the slow modes have the similar anisotropies and 
scalings. More promising are fast modes, which are isotropic (Cho & Lazarian 2002a). However, 
fast modes are subject to collisionless damping if the wavelength is smaller than the proton mean 
free path or by viscous damping if the wavelength is larger than the mean free path. According 
to CL02, fast modes cascade over time scales Tfk = T/ c x Va/v^ = (k x k m i n )~ 1 / 2 x Va/V 2 , where 
Tfc = kvk is the eddy turn-over time, V is the turbulence velocity at the injection scale. 

Consider gyroresonance scattering in the presence of collisionless damping. The cutoff of fast 
modes corresponds to the scale where Tfk'jd — 1 and this defines the cutoff scale k~ l . Using 
the tensors given in Eq. (29) we obtain the corresponding Fokker-Planck coefficients for the CRs 
interacting with fast modes by integrating Eq.(31) from k m i n to k c (see Fig.(7b)). When k~ l is 
less than the results of integration for damped and undamped turbulence coincides. Since the 
damping increases with (3, the scattering frequency decreases with (3. 

Adopting the tensors given in Eq. (30), it is possible to calculate the scattering frequency of 
CRs in high f3 medium. For instance, for density n = 0.5cm -3 , temperature T = 8000K, magnetic 
field Bq = 1/uG, the mean free path is smaller than the resonant wavelength for the particles with 
energy larger than 0.1 GeV, therefore collisional damping rather than Landau damping should be 
taken into account. Nevertheless, our results show that the fast modes still dominate the CRs’ 
scattering in spite of the viscous damping. 

Apart from the gyroresonance, fast modes potentially can scatter CRs by transit-time damping 
(TTD) (Schlickeiser & Miller 1998). TTD happens due to the resonant interaction with parallel 
magnetic mirror force —(mu^/2R)V||B. For small amplitude waves, particles should be in phase 
with the wave so as to have a secular interaction with wave. This gives the Cherenkov resonant 
condition u — knV\i ~ 0, corresponding to the n = 0 term in Eq.(31). From the condition, we see 
that the contribution is mostly from nearly perpendicular propagating waves (cos 6 ~ 0). According 
to Eq. (29),we see that the corresponding correlation tensor for the magnetic fluctuations My are 
very small, so the contribution from TTD to scattering is not important. 

Self-confinement due to the streaming instability has been discussed by different authors (see 
Cesarsky 1980, Longair 1997) as an effective alternative to scatter CRs and essential for CR accel¬ 
eration by shocks. However, we will discuss in our next paper that in the presence of the turbulence 
the streaming instability will be partially suppressed owing to the nonlinear interaction with the 
background turbulence. 

Thus the gyroresonance with the fast modes is the principle mechanism for scattering cosmic 
rays. This demands a substantial revision of cosmic ray acceleration/propagation theories, and 
many related problems may need to be revisited. For instance, our results may be relevant to the 
problems of the Boron to Carbon abundances ratio. We shall discuss the implications of the new 
emerging picture elsewhere. 
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6. Summary 

Recently there have been significant advances in the field of compressible MHD turbulence and 
its implications to cosmic ray transport: 

1. Simulations of compressible MHD turbulence show that there is a weak coupling between 
Alfven waves and compressible MHD waves and that the Alfven modes follow the Goldreich- 
Sridhar scaling. Fast modes, however, decouple and exhibit isotropy. 

2. Scattering of cosmic rays by Alfvenic modes is suppressed and therefore the scattering by fast 
modes is the dominant process provided that turbulent energy is injected at large scales. 

3. The scattering frequency by fast modes depends on collisionless damping or viscous damping 
and therefore on plasma (3. 
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Fig. 7— The scattering frequency v vs. the kinetic energy of cosmic rays (a) by Alfvenic 
turbulence, (b) by fast modes. In (a), the dash-dot line refers to the scattering frequency for 
isotropic turbulence. The ’x’ represents our numerical result for anisotropic turbulence, the solid 
line is our analytical result from Eq. (32). Also plotted (dashed line) is the previous result for 
anisotropic turbulence in Chandran (2001). In (b), the dashed line represents the scattering by 
fast modes not subjected to damping, the solid and dashdot line are the results taking into account 
collisionless damping. (From Yan & Lazarian 2002) 











